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ABSTRACT 
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We present optical spectra of the peculiar Type la supernova (SN la) 1999ac. The 
data extend from —15 to +42 days with respect to B-band maximum and reveal an event 
that is unusual in several respects. Prior to B-band maximum, the spectra resemble those 
of SN 1999aa, a slowly declining event, but possess stronger Si II and Ca II signatures 
(more characteristic of a spectroscopically normal SN). Spectra after B-band maximum 
appear more normal. The expansion velocities inferred from the Iron lines appear to 
be lower than average; whereas, the expansion velocity inferred from Calcium H and 
K are higher than average. The expansion velocities inferred from Si II are among the 
slowest ever observed, though SN 1999ac is not particularly dim. The analysis of the 
parameters vio(Si II), 7£(Si II), v, and Amis further underlines the unique characteristics 
of SN 1999ac. We find convincing evidence of C II A6580 in the day —15 spectrum with 
ejection velocity v > 16, 000 km s _1 , but this signature disappears by day —9. This rapid 
evolution at early times highlights the importance of extremely early-time spectroscopy. 

Subject headings: supernovae: general - supernovae: individual (SN 1999ac) 

1. Introduction 

Our understanding on the nature of Type la Supernovae (SNe la) would be greatly enhanced 
with the availability of high-quality spectra at very early phases (the first week or so after the 
explosion). At early phases, the ejecta is denser and observations probe the faster moving outer 
layers of the explosion (Ruiz-Lapuente et al. 1992; Khokhlov 1991). The diversity in Type la 
Supernova spectra is greatest at this time. At later phases, the ejecta thins, we probe deeper into 
the slower moving layers of the explosion and SN la spectra become more homogeneous. 

Mapping the composition of SN la ejecta at high velocity leverages direct constraints placed 
on explosion models. Such constraints include identifying the presence of ion signatures in SN 
spectra at various velocity intervals, a nontrivial task since the corresponding line profiles are often 
heavily blended. In particular, the presence of carbon at high velocity is consistent with a one- 
dimensional deflagration model (Nomoto et al. 1984; Branch et al. 1985). Silicon-peak elements at 
high velocity are more consistent with, for example, a one-dimensional delayed-detonation model 
(Khokhlov 1991). The signature of unburned material at a variety of velocities discovered across a 
sample of SNe la may support more recent three-dimensional deflagration models (Reinecke et al. 
2002; Gamezo et al. 2003). 

However, the task of obtaining early phase spectra presents numerous logistical challenges. 
Efficient discovery of SNe very shortly after explosion requires both high sensitivity and reliable 
rejection of false positives. A spectrum for confirmation (that a candidate is a SN) and classification 
must be obtained within a day. A program for intensive followup may then be triggered, requiring 
coordination of limited telescope time at possibly several different sites. Seeking to augment the 
local SN la dataset used as calibrated candles for cosmological distance measurement (Perlmutter 
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et al. 1998; Garnavich et al. 1998; Schmidt et al. 1998; Riess et al. 1998; Perlmutter et al. 1999; 
Knop et al. 2003; Tonry et al. 2003; Riess et al. 2004), the Supernova Cosmology Project (SCP) 
organized such a program together with several other teams (Aldering 2000) which ran during the 
Spring of 1999. When SN 1999ac was reported (Modjaz et al. 1999), the SCP had many telescope 
nights pre-scheduled for spectroscopy and photometry to observe this bright event. As a result, 
some of the earliest observed spectra of a SN la were obtained. 

SN 1999ac (R.A. = 16 h 07 m 15.0 s Decl. = +07 d 58'20", equinox 2000.0) was discovered and 
confirmed on unfiltered observations taken on Feb. 26.5 and 27.5 UT at 23". 9 east and 29". 8 south 
of the nucleus of its Scd host galaxy, NGC 6063 (Modjaz et al. 1999). Phillips et al. (1999a) reports 
that a confirmation spectrum taken on Feb. 28 UT is similar to that SN 1999aa with stronger 
Si II A6355 and well defined Ca II H&K. Optical light curves of SN 1999ac have been discussed 
in Phillips et al. (2002) and Li et al. (2003). They report a B-band light curve similar to that of 
SN 2002cx raising as fast as SN 1991T and declining similarly to SN 1994D until two weeks after 
maximum. The V-band light curve of SN 1999ac closely resembles that of SN 2002cx until 30 days 
after maximum light. The similarities between SN 1999ac and SN 2002cx break down in R and I 
bands. 

SN 1999ac was found to be spectroscopically similar to SN 1999aa (Li et al. 2001b; Garavini 
et al. 2004), which is considered to be an "intermediate" SN la with properties "between" those of 
spectroscopically normal SNe la and those of the spectroscopically peculiar, bright SN 1991T (for 
a discussion of whether or not this has implications for progenitor channels, consult Li et al. 2001b; 
Branch 2001). However, as we will show in this paper, the peculiarities of SN 1999ac are not limited 
to early epochs. The study of the ejecta geometry indicates calcium lines being formed in high 
velocity layers and iron and silicon lines being formed in low velocity layers. Furthermore, when 
SN 1999ac is plotted in the parameter space of 7£(Si II), «io(Si II), v, and Amis it stands out as an 
outlier among other well studied SNe. 

Here we present the collected spectra (our photometry of this SN will appear in a later work) 
and discuss in detail the two earliest spectra taken at days —15 and —9 with respect to the date 
of maximum brightness. The organization of this article is as follows. In §2, a description of the 
reduction scheme is given. Section §3 presents the spectra of SN 1999ac, and includes some empirical 
analysis of the early time spectra. Fits to the early spectra produced with the SN spectrum synthesis 
code SYNO W appear in §4. A study of the ejecta geometry is carried out in §5 and a comparison 
of several spectral indicators with those of other objects in §6 and §7 concludes the article. 

2. Data & Reduction 

Figure 1 shows the position of SN 1999ac in its host galaxy NGC 6063, an Scd galaxy with a 
recession velocity of 2848 km s _1 as determined from narrow Ha emission (Theureau et al. 1998). 
All spectra presented in this paper have been shifted to rest frame using this recession velocity. 
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The data set consists of 14 optical spectra extending from day —15 to day +42 (in this work all 
phases are expressed with respect to B-band maximum). In most cases, the spectra were acquired 
using different instrumental settings for the blue and red parts of the spectrum to avoid possible 
second-order contamination. Hence, the fully reduced 14 spectra are the combinations of both blue 
and red parts. The observation log appears in Table 1. 

The data were reduced using standard IRAF routines. The two-dimensional images were bias- 
subtracted and flat-fielded. The sky background was fitted, subtracted, and extracted for system- 
atics checks on the wavelength calibration. Wavelength and flux calibration were applied to the 
one-dimensional extracted spectra using calibration observations taken with the same instrumental 
setting and during the same night as science observations. The accuracy of the wavelength cali- 
bration was checked against the extracted sky spectra and generally found to agree to within 2 A. 
Synthetic spectrophotometry was computed to check the overall accuracy of the flux calibration 
against the publicly available photometric data (Jha 2002) and found to agree within the quoted 
statistical uncertainties. An atmospheric extinction correction was applied via tabulated extinction 
coefficients for each telescope used. 

Galactic reddening in the direction of SN 1999ac was estimated following Schlegel et al. (1998). 
A value of E(B — V) = 0.046 mag was found. For host galaxy extinction, Phillips et al. (2002) 
found that SN 1999ac's photometric evolution does not follow the Lira-Phillips relation (Phillips 
et al. 1999b), making problematic an estimate of host galaxy extinction with this technique. Labbe 
(2001) tentatively derive a total extinction to SN 1999ac of Ay = 0.51 mag. All spectra of SN 1999ac 
were corrected for Galactic extinction assuming Ry = 3.1, i.e. Ay = 0.14 (Cardelli et al. 1989) 
but no attempt was made to correct for host galaxy extinction correction because of the mentioned 
uncertainties in its estimate. This has no impact on the analysis presented, since we focus our 
analysis on line profiles, which cover wavelength scales smaller than those affected by reddening. 
For a given line profile, the effect of reddening can be considered as a global flux offset. 

The amount of host galaxy light contamination was estimated by x 2 -^, of the data with a SN 
spectrum template. It was found to be negligible at all epochs, and hence no subtraction of host 
galaxy light beyond that included in the sky background reference regions was performed. 

Telluric corrections are applied to spectra used for synthetic spectrum analysis, using calibration 
spectra from spectrophotometric standard stars in the vicinity of SN 1999ac. The wavelength regions 
with telluric corrections are marked in the figures, and their effect can be seen by comparing Figure 
2 with Figure 4, which show the spectra after telluric correction. 

3. Spectra 

The timespan and sampling frequency of the data set permits study of the spectroscopic evo- 
lution of SN 1999ac from very soon after explosion to seven weeks beyond B-band maximum light. 
Figure 2 presents the 14 fully reduced spectra, along with phases of observation for each. The top 
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spectrum, taken at day —15, is one of the earliest ever obtained of a SN la. 

The two earliest spectra are relatively featureless, though they clearly show the ion signatures 
of a SN la. The more obvious absorption features are at 4900, 6100, and 8100 A; due to Fe III 
blends, the Si II A6355 blend, and the Ca II infrared (IR) triplet (respectively). A weak absorption 
at 4200 A is probably also due to Fe III, while Si III is responsible for a weak absorption at 4400 A. 
Two very weak notches at 5200 and 5500 A hint at the presence of the S II "W" feature. A small 
depression just redward of the Si II 6100 A absorption could be due to C II. We return to this issue 
later in this article. 

The absorptions strengthen up to about day +8. The Si II A6355 absorption shifts noticeably 
to the red during this time. After day +8, the Ca II IR triplet profile morphology remains roughly 
constant. The S II "W" feature strengthens at maximum light, but shortly afterward disappears as 
it is replaced by Fe II and Na I lines. 

The absorption at 6100 A, due to Si II A6355, changes shape between days +8 and +24 as 
Fe II lines strengthen and obliterate it through line blending. By day +24, four robust minima have 
replaced the 6100 A absorption. Furthermore, the ostensible Si II emission peak appears to shift 
redward during this interval. The simplest explanation for this behavior is not a real shift per se, 
but rather an effect of blending, as the bluest part of the emission feature is overcome by the reddest 
absorption notch. It is interesting to note that Fe II lines begin contaminating the Si II feature as 
early as day +11, somewhat earlier than usual in SNe la. 

From day +24 onward, however, the spectrum basically does not evolve. The transition to 
iron-peak species dominance is complete just three weeks after maximum light, as numerous iron- 
peak lines dominate the spectrum from the ultraviolet to the near IR, excluding the Ca II H&K 
and IR triplet features. 

3.1. Early-Time Comparisons 

Generally, the post-maximum spectroscopic evolution of SN 1999ac resembles that of a normal 
SN la. The early-time spectra (days —15 and —9) present the unusual opportunity of comparison 
with the small set of supernovae with similarly early spectra. 

In Figure 3, the day —15 spectrum is presented along with early-time spectra of SNe 1991T 
(Filippenko et al. 1992), 1999aa (Garavini et al. 2004), 1990N (Leibundgut et al. 1991), and 1994D 
(Patat et al. 1996). Compared to the spectroscopically peculiar SNe 1999aa and 1991T, SN 1999ac 
exhibits stronger Ca II and Si II A6355 absorptions. On the other hand, its Fe III 4400, 4100 A and 
Si III A4560 features are weaker. 

A small flux depression in the day —15 spectrum is clearly visible to the red of the Si II 6100 A 
feature, echoing one visible in the spectrum of SN 1990N. Though made uncertain by the presence 
of the telluric absorption at 6900 A, another feature at 7000 A could be common to both spectra as 
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well. The first feature has been identified as C II A6580 in SN 1990N (Mazzali 2001), in SN 1994D 
(Hatano et al. 1999a), in SN 1998aq (Branch et al. 2003), and in SN 1999aa (Garavini et al. 2004). 
The second feature could be due to another line from the same ion, C II A7234. 

Overall, the —15 day spectrum resembles that of SN 1990N more than it does SN 1994D, with 
its more obvious (and characteristic) S II and Fe II blends. In summary, at day —15, SN 1999ac 
possess spectroscopic characteristics common to both slow-declining SNe like SN 1991T and spec- 
troscopically normal ones like SN 1990N at a similar phase. 

The —9 day spectrum of SN 1999ac (Figure 3) does not extend to the Ca II H&K absorption. 
The two Fe III features are still weaker than in SN 1991T and possess rounded minima. The 
contribution from Fe II in the region between 4000 and 5000 A appears weaker than in SN 1994D 
and in SN 1990N. The absorption feature present on the red edge of this line and that around 
7000 A at day —15 are no longer evident. 

4. Synthetic Spectra 

One of the more interesting ion signatures suggested in the previous section is that of C II. The 
presence of carbon lines and their ejection velocities pose potential constraints on hydrodynamical 
explosion models. Carbon at high velocity is consistent with both one-dimensional (Nomoto et al. 
1984) and three-dimensional (Reinecke et al. 2002; Gamezo et al. 2003) deflagration models. An 
apparently generic result from the newer models is the mixing of "fuel" (carbon and oxygen) and 
"ashes" (products of nucleosynthesis) at all velocities. Hence, the detection of carbon at low velocity 
would support the three-dimensional deflagrations. 

To further explore the C II signature and its behavior at early times, in this section we compare 
synthetic spectra to those observed before maximum light. Clearly, the ideal approach would be to 
invert a SN spectrum in some way to yield a composition model. Unfortunately, the SN atmosphere 
problem is an ill-posed inverse problem. Instead, one of two approaches must be adopted. 

One approach is detailed analysis, where the goal is to include all relevant transfer physics 
(nonlocal thermodynamic equilibrium rates, relativity, time-dependence, energy from radioactive 
decays, etc.) and numerically simulate the emergent spectrum of a given hydrodynamical model. 
This approach is appropriate for validating hydrodynamical models, or suggesting adjustments to 
such models in the future. Though powerful, detailed analysis codes (for example, the general- 
purpose PHOENIX code, Hauschildt & Baron 1999) consume months of computer time for a single 
calculation. 

Another, complimentary approach is direct analysis. The goal of direct analysis is more empir- 
ical; to constrain the presence or absence of ions in a spectrum and the ejection velocities of their 
parent atoms. The task is generally nontrivial, since SN spectra consist of many blends of lines 
which cannot be treated simplistically. Still, approximate techniques are used to make the process 
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fast and iterative, but the results are powerful. Constraints from direct analysis are of use to both 
detailed modelers and explosion modelers. In the latter case, the results of direct analysis can rule 
out many hydrodynamical models (in principle) before the extensive simulation calculations are 
performed. 

SYNOW (Fisher 2000) is a direct analysis code that generates spectra based on a simple, 
conceptual model of a SN appropriate during the first few weeks to months after explosion. This 
model consists of a blackbody-emitting, sharply defined photosphere surrounded by an extended line- 
forming, pure scattering atmosphere. The entire envelope is assumed to be homologously expanding. 
Line transfer is treated using the Sobolev method (Sobolev 1960; Castor 1970; Jeffery & Branch 
1990) so line opacity is parameterized in terms of Sobolev optical depth. Which ions are used in 
the calculation is determined by experience, guided by the SN ion signatures atlas of Hatano et al. 
(1999b). For each ion introduced, Sobolev optical depth as a function of velocity for a "reference 
line" (usually a strong optical line) is specified. Optical depths in other lines of the ion are set 
assuming Boltzmann excitation of the levels at temperature T exc . 

The parameters v p hot and set the velocity and blackbody continuum temperature of the 
photosphere, respectively. For each ion, the optical depth r at the minimum ejection velocity v m i n 
is specified. Optical depth scales exponentially with velocity according to the e-folding velocity 
parameter v e , and is considered to be zero for velocity greater than v max . If v m i n > v p hot for an ion, 
we refer to the ion as "detached." 

A sharply defined, blackbody-emitting photosphere clearly cannot serve as a perfect substitute 
for the processes of continuum formation in a SN atmosphere. Hence, the synthesized continuum 
level may systematically differ from that observed in some wavelength regions. Generally, a good 
fit to the blue continuum results in a brighter red synthetic continuum. This effect has little to no 
bearing on line identifications or velocity inferences. 



Day -15 

Figure 4 is a comparison of a synthetic spectrum to the observed spectrum at day —15. Choosing 
Tbb = 11, 200 K reproduces a satisfactory reproduction of the overall continuum shape. We also find 
that the choice of v p hot = 13, 000 km s _1 assists in producing reasonable line profiles when other 
factors are taken into account for each individual ion. 

The Si II optical depth profiles are detached, with v m i n = 14,200 km s _1 . This detachment 
moves the synthetic Si II A6355 absorption higher in velocity space so that it matches the observed 
feature. It also flattens the corresponding emission feature to improve the agreement. We regard 
the presence of Si II in this spectrum as definite. 

Si III optical depth must be capped at v max = 17,000 km s _1 to match the feature usually 
associated with it at 4400 A. Interestingly, a concomitant synthetic feature appears to reproduce 



9 



a feature in the near-infrared around 9000 A. Though the fit is not perfect here, we regard the 
presence of Si III in this spectrum as definite. 

The velocity range in which we introduce Ca II is mainly constrained by the Ca II IR triplet 
since the H&K component is missing from the data. Nevertheless, the agreement with the observed 
Ca II IR triplet and the falling edge of the observable H&K signature is convincing. The presence 
of Ca II in this spectrum is definite. 

Introducing Fe III provides a match to two observed features, one at 4200 A and another 
at 4900 A. Like Si II, the corresponding optical depth profile must be detached (in this case to 
Vmin = 14,500 km s _1 ) to reproduce the features. Adding some Mg II improves the fit to the 
absorption at 4200 A. The presence of Fe III is definite. 

As previously mentioned, the characteristic S II feature does not yet appear to be fully developed 
at phase —15 days. Including some S II optical depth produces an absorption blend around 5300 
A. Generally, however, fits to the S II blend region using SYNOW do not reproduce the observed 
features in most cases. Our fit in this region is not optimal, nevertheless, we regard the presence of 
S II in the day —15 spectrum as probable. 

It is interesting that including some Ni III to the synthetic spectrum helps to account for the 
flux deficit to the blue of the Fe III 4900 A feature and near 5300 A. The appearance of this ion in a 
spectrum is somewhat unusual, nevertheless the improvement cannot be discounted. Conservately, 
however, we regard the presence of Ni III as possible. 

The weak absorption at 6300 A is well matched by a detached C II (v m i n = 16,000 km s _1 ) 
that may also contribute to a feature near 4500A and perhaps at 7000 A. The wavelength regions 
where C II makes its contribution are highlighted in Figure 5 showing the effect of the presence of 
this ion on the synthetic spectrum. The good matching of the absorption feature at 6300 A and the 
possible contribution near 4500 A makes the identification of C II definite in the day —15 spectrum. 

We have considered C III to match the small notch on the red side of the Si III feature near 
4500 A. However, we find the evidence for C III less convincing than that for C II. The small optical 
depth used to generate the C III feature prevents the appearance of other, weaker C III lines in the 
synthetic spectrum. Hence, the only sign of C III (an unusual identification) would be this line by 
itself. Without concomitant C III lines, we consider this ion to be only a remote possibility. The 
feature is probably not due to H/3, since the spectrum shows no other signs of hydrogen. 

Day -9 

Figure 6 presents the synthetic spectrum produced for comparison with the day —9 spec- 
trum. The fit parameters used appear in Table 3. The photospheric velocity has been lowered to 
Vph ot = 11, 800 km s _1 , and the blackbody temperature increased to = 13, 800 K. The increased 
temperature should not be interpreted as a peculiarity since the blackbody continuum used by 
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SYNOW has to be regarded only as an approximation with little physical meaning (see section 4 
for details.) For the most part, the same ions are used for this fit as in the previous one. Note 
that at this phase, SYNOW is overestimating the continuum through the Si II 6100 A absorption 
and redward. However, this offset has no effect on line identifications or ejection velocity interval 
measurements. 

Ions definitely present in the spectrum include Si II (though no longer detached), Si III and 
Fe III. Again the fit to S II is problematic, but its identification based on the two notches redward 
of 5200 A is not unreasonable. Ni III seems to remove some excess flux again near 4700 A, but it 
no longer is sufficient to account for all of the flux deficit. 

Most interestingly, it appears that between day —15 and day —9, the C II signature has 
disappeared. A small notch near 4500 A could be due to a small amount of C III at photospheric 
velocities, but adding C III optical depth only enhances the excess flux to the red of this feature. 
For this reason, we consider the presence of C III in this spectrum as unlikely, though no alternative 
has been identified. 



5. Ejecta geometry 

We have seen in previous sections tentative evidence of C II lines moving at velocities above 
20,000 km/sec. The peculiarity of SN 1999ac is not limited to the high velocity ejecta. Distinctive 
characteristics are noticed in many other absorption features pointing out an overall unusual ejecta 
geometry. This can be investigated, in first approximation, by looking at the line profile of spectral 
features. The characteristic line profile during the first several weeks after outburst is a P-Cygni 
profile. We here briefly review the relevant features of the P-Cygni line profile before presenting 
measurements based on it. 

A P-Cygni profile arises from a configuration consisting of an extended, expanding line-forming 
region surrounding an optically thick core. Consider such a configuration, observed in the frame 
coinciding with the rest frame of the core's center. Material in front of the core (as observed) 
moves toward the observer, and scatters radiation out of the line of sight, resulting in an absorption 
feature blueshifted with respect to the line rest wavelength. Material not in front of the core (and 
expanding away from it) scatters radiation into the observer line of sight, resulting in an emission 
feature peaked at the line rest wavelength. Together, the blueshifted absorption and rest-wavelength 
centered emission features form the P-Cygni profile. 

The morphology of a P-Cygni absorption component provides an estimate of the velocity inter- 
val in which the originating line forms. The blue edge of the absorption feature forms at the highest 
ejection velocities where the line is optically thick. In practice, the minimum of an absorption profile 
is used to derive a characteristic ejection velocity describing where in velocity space a line forms. 

The strength of a line can easily influence its shape. A weak line (with smaller line opacity) 



- 11 - 



produces a sharp, robust minimum that can be measured with little ambiguity. Conversely, a strong 
line produces a more rounded absorption feature, and as line strength increases, the position of the 
absorption minimum shifts to the blue. Hence, a velocity measured from such a minimum is not 
necessarily representative of the minimum ejection velocity of a line. 

Two lines from the same parent ion, or concomitant lines, that are closely spaced in wavelength 
(as in a doublet) may generally be treated as a single line for the purposes of velocity estimation. On 
the other hand, features from other ions blending with a given line present a problem, as they can 
easily shift the position of the line minimum to the blue or red. Therefore, care must be exercised 
when making inferences from absorption features that may or may not be suffering blending effects. 

5.1. Ejecta geometry comparison 

In the following paragraphs we compare the velocity field of some of the characteristic features of 
SN 1999ac with those of other SNe known for having a peculiar ejecta geometry, namely SN 2002cx 
(Li et al. 2003), SN 1999aa (Garavini et al. 2004) and SN 2002bo (Benetti et al. 2004b), and with that 
of SN 1994D (Patat et al. 1996) which we here take as a prototype of normal supernova. SN 2002cx 
is a well studied under-luminous supernova with normal B-band light curve decline rate, spectral 
signatures similar to SN 1991T but with spectral lines with low expansion velocities. SN 2002bo 
had a normal decline rate but showed spectra with higher than average expansion velocities similar 
to SN 1984A (Barbon et al. 1989). Finally, SN 1999aa was found to have a slow light curve decline 
rate and very weak Si II absorption features with expansion velocities constant in time. Also, its 
spectra had similarities with SN 1991T at early epochs, but rapidly changed toward normal looking 
spectra just before maximum light. 

5.1.1. P re-maximum spectra 

The comparison of the spectrum of SN 1999ac at day —9 with those of the SNe mentioned 
above is shown in Figure 7, Panel A. The most noticeable differences lie in the region of Fe III and 
Si II lines. To investigate such differences, we show, in Panel B and Panel C, the comparison in 
velocity space of the absorption features in the wavelength region respectively around 4250 A and 
6150 A. 

In early spectra the absorption visible at 4250 A is generally due mostly to Fe III A4404 but 
a contribution of Mg II A4481 is expected and can vary from object to object. Normal SNe have 
a stronger Mg II component whereas SN 1991T-like SNe are dominated by Fe III. Among the SNe 
we compare, the Mg II contamination is probably more important for SN 1994D and SN 2002bo. 
SN 1999ac appears to have approximately the same velocity distribution as SN 2002bo and as 
SN 1999aa. Both SN 1994D and SN 2002cx have lower velocities than our object. The actual ejecta 
geometry is difficult to disentangle from this comparison because of the uneven contamination 
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of Mg II among the different objects. However, based on the analysis we performed in §4, this 
absorption in SN 1999ac appears to be the result of both Mg II and a dominant component of low 
velocity Fe III. We will come back later to this point for more discussion. 

In Panel C we compare the line profile in velocity space of the absorption feature at around 
6150 A. This line is due to Si II A6355 and it is usually the one suffering the least blending among 
all the supernova optical spectral features. It is thus considered the simplest to use for expansion 
velocity studies. At this epoch SN 1999ac appears to have its line minimum at the same velocity 
as SN 1994D but with a lower velocity for its blue edge. This is consistent with the optical depth 
of Si II being weak in the outermost layer. We will come back to this point later in the analysis. 
SN 2002bo and SN 2002cx show respectively a faster and slower Si II layer. 

5.1.2. Spectra at maximum 

The ejecta geometry of the same SNe at around maximum light is compared in Figure 8, Panel 
A . As in the previous epoch the most noticeable differences are in the iron and silicon absorption 
features. Furthermore, Ca II H&K also appears to have a different line profile in each supernova. 

The line profile in velocity space of Ca II H&K (Panel B) of SN 1999ac is comparable with that 
of SN 2002bo and has a single minimum while SN 1994D and SN 1999aa have double minima ( see 
Nugent et al. (1997); Garavini et al. (2004) a discussion of the double minima in these SNe.). 

At this epoch the velocity of the minima are approximately all comparable with differences 
of the order of a thousand kilometer per second. The only exception is SN 2002cx for which this 
feature is much weaker and at lower velocity. 

Fe II A5083.4 (Panel C) is at lower velocity in SN 1999ac than in SN 1994D, SN 2002bo and 
SN 1999aa, while only SN 2002cx shows even lower values. This is consistent with the iron layer of 
SN 1999ac being deeper into the atmosphere than that of the other objects as already mentioned 
for the case of Fe III A4404. At this epoch Si II A6355 (Panel D) also shows lower values than all 
the other SNe analyzed except for SN 2002cx where this line is considerably weaker. 

5.1.3. Late time spectra 

In late time spectra, showed in Panel A of Figure 9, the absorption features are mainly formed 
by iron lines. The velocity distribution and relative abundances of iron in the SNe analyzed can be 
discerned by inspection of the different strengths of the small absorptions and peaks in the spectra. 
In Panel C, vertical dotted lines mark the position of the various iron line minima. SN 1999ac 
shows lower velocity compared to all SNe with the usual exception of SN 2002cx. Si II at this epoch 
is too weak and blended with iron lines to disentangle its complete line profile and is not analyzed 
further. Ca II H&K (Panel B) is still strong, and for SN 1999ac this absorption feature has the 
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highest velocity among the analyzed SNe; The Ca II IR, Panel D, has a velocity similar to the other 
analyzed SNe. 

5.2. Velocity time evolution 

The most striking peculiarity of SN 1999ac seems to be the concomitance of high velocity 
calcium with low velocity iron and silicon. This characteristic becomes more evident with time. 
Table 4 lists the velocities inferred from the minimum of the Ca II H&K feature, together with 
those measured for SN 1999aa (Garavini et al. 2004). Figure 10 shows these velocities as a function 
of time compared with those of other well observed SNe (see Garavini et al. (2004) for details on the 
measurement techniques). The statistical uncertainty of the measurements is negligible compared 
with the estimated 100 km s _1 uncertainty in host galaxy recession velocity. Spectra obtained before 
maximum light do not include this feature, so at those phases no measurement can be reported. 
After maximum, the Ca II H&K velocity decreases monotonically. This is in agreement with the 
trend measured for most of other SNe, but with SN 1999ac systematically offset to higher velocity. 

The evolution of the velocity inferred from the minimum of Si II A6355 in Figure 11 is even 
more peculiar. (The measured values are reported in Table 4 together with those measured for 
SN 1999aa.) The Si II velocity of SN 1999ac appears to be monotonically decreasing with time 
as is usually the case for dimmer supernovae, like SN 1991bg or SN 1999by. However, SN 1999ac 
was not intrinsically dim (Mg a:E =-18. 98(39) Li et al. (2003)). Some events, such as SN 2000cx or 
SN 1999aa, maintain a rather constant Si II velocity as a function of time. It should be noted that 
after day +11 the measurements may become less meaningful because of possible blending with 
Fe II. Figure 12 illustrates this ambiguity clearly: By day +24, the entire feature can no longer be 
assumed to arise from Si II A6355 alone, since the peak of the emission feature has shifted to the 
red of the rest wavelength (zero velocity). Still, it is clear that the velocity evolution of the feature 
out to day +11 is extreme among SNe la. 

6. Type la SN parameter space. 

Type la supernovae are currently believed to be a multiparameter class of objects. The standard 
paradigm is, however, to describe the intrinsic spread of properties through one parameter: the light 
curve width stretch, s, or equivalently, through the light curve decline rate Amis. This, indeed has 
proved to be successful in reducing the intrinsic spread in brightness and in enabling the use of Type 
la SNe as high-quality distance indicators. Several other parameters have been proposed with the 
goal of fully describing SNe la and their parameter space. In the following sections we compare the 
measurements of the parameters 7£(Si II) (Nugent et al. 1995), i>io(Si II) (Branch & van den Bergh 
1993) and v (Benetti et al. 2004a) for SN 1999ac with those of the dataset presented in Benetti 
et al. (2004a). For completeness we report here also the data for SN 1999aa analyzed in Garavini 
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et al. (2004). The measured values for both supernovae are reported in Table 5. The symbols used 
in the following figures are chosen to be as those in Benetti et al. (2004a) to maintain the same 
distinction among the three clusters found in their analysis. 

6.1. Ami 5 versus K(Si II) 

Figure 13 is a plot of Amis as a function of the quantity TZ(Si II) for many SNe la, including 
SN 1999ac. Nugent et al. (1995) defined fc(Si II) as the ratio of the depth of the Si II 5800 A 
absorption to that of the Si II 6100 A absorption. The authors theorize that the observed correlation 
between Amis and 7£(Si II) is driven by temperature (and hence nickel mass). Thus, hotter and 
brighter events tend to be characterized by a small TZ(Si II) value; cooler and dimmer events are 
characterized by a larger 7£(Si II) value. The plot shows that SN 1999ac is above the general 
trend. The small value of 7£(Si II) of SN 1999ac is consistent with the presence of strong Fe III 
lines in early epochs and a high temperature of the envelope. The photometric peculiarities of 
SN 1999ac, reported in section 1, suggest that for this object Ami 5 might not be a good indicator 
of luminosity. The parameter space of this diagram is further expanded by the inclusion of the 
values for SN 1999aw, SN 1999bp and SN 1999aa (Folatelli 2004) to which the correlation appears 
to apply. 

6.2. Amis versus t>io(Si II) 

Hatano et al. (2000) investigated the spectroscopic diversity of SNe la by plotting values of 
7£(Si II) against vio(Si II) (the blueshift measured in the Si II A6355 feature at ten days after 
maximum light). The authors reasoned that if SNe la were a one-dimensional family based on the 
mass of synthesized 56 Ni, then the two observables would be correlated. Unable to discern such a 
correlation, they reasoned that the differences reflect differences in the explosion mechanism itself. 
In Figure 14, we plot Amis (a proxy for 7£(Si II)) against t>io(Si II), including observed values 
from SN 1999ac, SN 1999aw, SN 1999bp and SN 1999aa. SN 1999ac possesses a low value for 
vio(Si II) but a normal value of Amis (Li et al. 2003), similar to that of SNe 1989B, 1994D, and 
1996X, making SN 1999ac unique in falling outside the apparent groupings. The same would be 
true plotting 7£(Si II) instead of its proxy Amis. While SN 1999ac is spectroscopically similar to 
SNe 1990N and 1999aa (except with lower velocities), its placement on the diagram away from those 
events indicates the possibility of even higher dimensions of diversity in these objects. 

The parameter space of this diagram is further expanded by the inclusion of the values for 
SN 1999aw, SN 1999bp and SN 1999aa, though perhaps in directions that appear more consistent 
with the trends among the groupings. 
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6.3. Amis versus v 

The parameter v was introduced in Benetti et al. (2004a) as an estimate of the expansion veloc- 
ity time derivative computed after B-band maximum light. The authors found a weak correlation 
with Amis, Fast declining, under- luminous supernovae show large v. Slow declining supernovae 
show small v, while normal supernovae can have both large and small values. In the Amis versus 
v plane (Figure 15) SN 1999ac falls on the high edge of the normal supernovae with the highest v 
measured to date - similar to that of SN 1983G, which had different t>io(Si II) and 7£(Si II). The 
parameter space of this diagram is again expanded by the inclusion of the values for SN 1999aa and 
SN 2000cx which are found to have very small values of v (Garavini et al. 2004). 

6.4. v versus vio(Si II) 

The plot of v versus t>io(Si II) shown in Figure 16, was left unexplored in Benetti et al. 
(2004a) but it is interesting to note that the three SN groups identified by the authors still populate 
different regions of this plot, meeting each other at intermediate values. On this plot underluminous 
supernova tend to have small values of v io(Si II) and high values of v. Normal supernovae populate 
the central part of the plane while supernovae with fast expansion velocity (i.e. high v io(Si II)) have 
also high v. SN 1999ac has the highest v and the lowest vio(Si II) among the supernova measured, 
making it similar to the faint supernovae but on the extreme end. However, as noted, SN 1999ac 
did not appear to be dim. It also should be noted that SN 2000cx falls in an unpopulated region of 
the plot with high uio(Si II) and low v. 

7. Conclusions 

We have presented spectroscopic observations of SN 1999ac from —15 to +42 days with respect 
to B-band maximum light. The earliest spectra are similar to those of SN 1999aa, but share some 
characteristics of spectroscopically normal SNe like SN 1990N. Notable is the early appearance of 
iron features in the spectrum at +11 days and the lack of any real evolution in the features after 
+24 days. 

Using synthetic spectra, we have unambiguously identified the presence of C II at high velocity 
in the earliest spectrum at —15 days. However, just six days later, all significant traces of C II 
have disappeared. This indicates that obtaining spectra even earlier than —10 days with respect 
to maximum may be required to reliably probe the outer layers of SNe la. Obtaining this level 
of efficiency from an observing program requires careful coordination, similar to that achieved by 
the European Research Training Network 1 search for SNe la, and the Nearby Supernova Factory 
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(Aldering et al. 2002). 

Comparing the spectra of SN 1999ac with those of other supernovae we find an unusual ejecta 
geometry. Iron and silicon lines appear to be formed in deeper atmosphere layers than the corre- 
sponding lines for other supernovae with the exception of the extreme case of SN 2002cx. Calcium 
lines, however, are found to be formed in high velocity layers. The same trends are confirmed when 
analyzing the time evolution of velocities as derived from the minimum of Ca II H&K and Si II 
A6355. The former shows a trend consistent with normal SNe la, though the values are slightly 
higher than average, while the latter shows monotonically decreasing values that follow the trend 
of under-luminous SNe la. 

We have presented measurements for spectroscopically derived observables such as 7£(Si II), 
uio(Si Il)and v. While the values for SN 1999ac weakly support the correlation between TZ(Si II) 
and Amis, the position of SN 1999ac on the Amis versus vio(Si II) plane is somewhat off the main 
trend. In the plane v versus Amis SN 1999ac falls in a region of normal supernovae but has the 
highest value of v, while in the plane v versus t>io(Si II) it appears to be on the extreme of faint 
supernovae with the highest v and lowest fio(Si II). However, SN 1999ac was not reported as a dim 
supernova. This analysis points out that SN 1999ac is unlike any other known supernova. 
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Fig. 1— SN 1999ac in its host galaxy NGC 6063. R.A. = le'W^lS.O 5 Decl. = +07 d 58 m 20 s 
(equinox 2000.0) 
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Fig. 2. — SN 1999ac spectral evolution. Epochs referred to B-band maximum light. The © symbol 
marks sharn atmosnheric absorntion features. 
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Fig. 3.— The -15 and -9 days spectra of SN 1999ac together with those of SN 1999aa, SN 1991T, 
SN 1990N and SN 1994D respectively from Garavini et al. (2004); Filippenko et al. (1992); Lei- 
bundgut et al. (1991); Patat et al. (1996). Epochs are quoted in the labels. Line identification are 
taken as in Li et al. (1999, 2001a); Fisher et al. (1999); Patat et al. (1996); Mazzali et al. (1995); 
Kirshner et al. (1993); Jeffery et al. (1992) 
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Fig. 4. — Synthetic spectra compared with SN 1999ac spectrum for —15 days. Dashed line: best 
match synthetic spectrum. Solid line: data. SYNOW parameters used are presented in table 2. 
Ions responsible for features in the synthetic spectrum are marked. Where a telluric feature has 
been removed, the spectrum is marked with an Earth symbol in parenthesis. 
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Fig. 5. — Synthetic spectra (heavy solid line) compared with SN 1999ac spectrum (light solid line) 
for —15 days in the 4500 A (left panel) and 6150 A (right panel) region. First model from the top: 
Ten = 0; Second model from the top: Ten / 0. SYNOW parameters used are presented in table 2. 
Ions responsible for features in the synthetic spectrum are marked. 
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Fig. 6. — Synthetic spectra compared with SN 1999ac spectrum for —9. days Dashed line: best 
match synthetic spectrum. Solid line: data. SYNOW parameters used are presented in table 3. 
Ions responsible for features in the synthetic spectrum are marked. Where a telluric feature has 
been removed, the spectrum is marked with an Earth symbol in parenthesis. 




Fig. 7.— Panel A: The -9 day spectrum of SN 1999ac together with those of SN 1999aa, SN 2002cx 
and SN 1994D. Epochs are quoted in the labels. Line identifications are taken as in Li et al. (1999, 
2001a); Fisher et al. (1999); Patat et al. (1996); Mazzali et al. (1995); Kirshner et al. (1993); Jeffery 
et al. (1992). Panel B: Comparison of the same SNe in the region of Fe III A4404 in velocity space. 
Panel C: Comparison of the same SNe in the region of Si II A6355 in velocity space. In Panel B 
and C vertical dotted lines indicate the position of the feature minima for SN 1999ac. 
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Fig. 8. — Panel A: The spectrum at maximum light of SN 1999ac together with those of SN 1999aa, 
SN 2002cx and SN 1994D. Epochs are quoted in the labels. Line identifications are taken as in Li 
et al. (1999, 2001a); Fisher et al. (1999); Patat et al. (1996); Mazzali et al. (1995); Kirshner et al. 
(1993); Jeffery et al. (1992) Panel B: Comparison of the same SNe in the region of Ca II H&K in 
velocity space. Panel C: Comparison of the same SNe in the region of Fe II A5083.4 in velocity 
space. Panel D: Comparison of the same SNe in the region of Si II A6355 in velocity space. In Panel 
B, C and D vertical dotted lines indicate the position of the feature minima for SN 1999ac. 
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v[km/s] Rest Frame X[A] v|km/s] 

Fig. 9.— Panel A: The +24 day spectrum of SN 1999ac together with those of SN 1999aa, SN 2002cx 
and SN 1994D. Epochs are quoted in the labels. Line identifications are taken as in Li et al. (1999, 
2001a); Fisher et al. (1999); Patat et al. (1996); Mazzali et al. (1995); Kirshner et al. (1993); Jeffery 
et al. (1992) Panel B: Comparison of the same SNe in the region of Ca II H&K in velocity space. 
Panel C: Comparison of the same SNe in the Fe II lines region in wavelength space. Panel D: 
Comparison of the same SNe in the region of Ca II IR triplet in velocity space. In Panel B, C and 
D vertical dotted lines indicate the position of the feature minima for SN 1999ac. 
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Fig. 10. — Expansion velocities as measured from the minima of Ca II H&K of SN 1999ac compared 
with the values of other SNe taken from Wells et al. (1994); Garnavich et al. (2004); Kirshner et al. 
(1993); Patat et al. (1996); Jha et al. (1999) and references therein. Values for SN 1999ac are 
marked as filled circles. Measured values are reported in Table 4 together with those measured for 
SN 1999aa. 
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Fig. 11. — Expansion velocities as measured from the minima of Si II A6355 of SN 1999ac compared 
with the values of other SNe taken from Li et al. (1999, 2001a); Garnavich et al. (2004); Salvo et al. 
(2001) and references therein. Values for SN 1999ac are marked as filled circles. Measured values 
are reported in Table 4 together with those measured for SN 1999aa. 
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Fig. 12. — Evolution of the Si II A6355 feature in Doppler space. The velocity at the minimum of 
the absorption feature is at 13,000 km s _1 at —15 days, and moves to 8,000 km s _1 at +11 days 
relative to maximum. By day +16, contamination by Fe II lines has clearly begun. 
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Fig. 13. — Light curve decline rate Amis versus 7£(Si II) for the SNe la in Benetti et al. (2004a). 
SN 1999ac generally supports the trend toward higher luminosity at lower 7£(Si II) but with a 
smaller 7£(Si II) with respect to normal SNe. Measured values are reported in Table 5 together with 
those mesured for SN 1999aa. 
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Fig. 14. — Plot of Amis versus uio(Si II). SN 1999ac falls in a relatively unpopulated region of 
the plot with low value of uio(Si II). Measured values are reported in Table 5 together with those 
measured for SN 1999aa. 
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Fig. 15. — Plot of Ami5 versus v. SN 1999ac generally falls in the region of the plot where other 
normal SNe lie but has the highest measured v. Measured values are reported in Table 5 together 
with those measured for SN 1999aa. 
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Fig. 16. — Plot of v versus t>io(Si II). SN 1999ac falls in a relatively unpopulated region of the plot 
with the highest v and lowest t>io(Si II). SN 2000cx falls in a unpopulated region of the plot with 
low v and high uio(Si II). Measured values are reported in Table 5 together with those measured 
for SN 1999aa. 
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Table 1: Data set specifications. 



JD 


Epoch f 


Telescope 


Instrument 


A Range g 


FWHM a ' 9 


(S/N) b 


Comments 


-2400000 


rsf B max 






[A] 


[A] 






51236.89 


-15 


APO 


DIS 


3703-10307 


n 


125 


5696.16 c 


51240.95 


-9 


MDM 2.4m 


MARK III 


3827-8860 


16 


63 


e 


51251.39 


+0 


ESO 3.6m 


EFOSC2 


3331-7495 


32 


325 


d 


51253.84 


+2 


CTIO 4m 


RCSP 


3235-9263 


5 


76 


d 


51253.72 


+2 


NOT 


ALFOSC 


3285-9655 


34 


176 


5852.06 c 


51259.88 


+8 


CTIO 4m 


RCSP 


3227-9254 


6 


94 


d 


51262.89 


+11 


CTIO 4m 


RCSP 


3254-9278 


6 


114 


d 


51267.87 


+16 


CTIO 4m 


RCSP 


3239-9241 


6 


90 


d 


51275.98 


+24 


KPNO 4m 


T2KB 


3029-10401 


11 


59 


d 


51279.85 


+28 


ESO 3.6m 


EFOSC2 


3341-10255 


32 


237 


7440.45 c ' d 


51282.90 


+31 


Lick 3m 


KAST 


3321-10483 


5 


46 


5489.40 c 


51284.84 


+33 


ESO 3.6m 


EFOSC2 


3392-10128 


32 


195 


7363.23 c ' d 


51290.84 


+39 


ESO 3.6m 


EFOSC2 


3344-10194 


32 


146 


7435.50 c ' d 


51293.97 


+42 


Lick 3m 


KAST 


3268-8002 


5 


49 


5417.33 c 



"Average over wavelength. 

'Average signal-to-noise ratio per wavelength bin. 

c Beginning of red channel, [A]. 

^Negligible 2 nd order contamination. 

e Possible 2 nd order contamination above 7500 A. 

•^Rest Frame days. 

9 Rest Frame wavelength. 
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Table 2: Synow parameters for —15 days. The fit is shown in Figure 4. v p h ot =13000 km s 
T bb = 11200 K. 



Ion 


r 


Vmin 

10 3 kms" 1 


Vmax 

10 3 kms" 1 


T 

- 1 exc 

10 3 K 


10 3 kms" 1 


C ii 


0.038 


16 


40 


15 


5 


C in 


0.2 




14.2 


15 


5 


i 


0.2 




40 


15 


5 


Mg ii 


0.15 




40 


15 


5 


Si ii 


0.65 


14.2 


40 


15 


5 


Si in 


0.42 




17 


15 


5 


S ii 


0.2 




17 


15 


5 


Ca ii 


1.5 


16 


40 


15 


5 


Fe in 


0.55 


14.5 


18 


12 


5 


Co ii 


0.006 




40 


15 


5 


Ni in 


5 




40 


12 


5 



- 38 - 



Table 3: Synow parameters for —9 days. The fit is shown in Figure 4.v p /i O t=11800 km s , = 
13800 K. 



Ion 


r 


Vmin 


"max 


T 

- 1 exc 








10 3 kms" 1 


10 3 kms- 1 


10 3 K 


10 3 kms- 1 


C II 


0.015 


16 


40 


12 


5 


C in 


0.75 




12.8 


12 


5 


i 


0.1 




40 


12 


5 


Mg II 


0.2 




40 


12 


5 


Si ii 


1.2 




40 


12 


5 


Si in 


0.6 




16 


12 


5 


S ii 


0.2 




17 


12 


5 


Ca ii 


1.5 


16 


40 


12 


5 


Fe in 


0.65 




18 


12 


5 


Ni III 


7.0 




40 


12 


5 


Co II 


0.045 




40 


12 


5 
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Table 4: Measurements of the expansion velocity inferred from Ca II H&K and Si II A6355 for 
SN 1999ac and SN 1999aa (Garavini et al. 2004). The uncertainties on the velocity measurements 
are assumed to be 100 km s _1 , for further information see §5.2. 



SN 1999ac 


SN 1999aa 


Epoch 


Ca II H&K 


Si II 


Epoch 


Ca II H&K 


Si II 


days 


[km s _1 ] 


[km s" 1 ] 


days 


[km s _1 ] 


[km s _1 ] 


-15 




13409 


-11 


20155 




-9 




12169 


-7 




10519 





14045 


9801 


-3 


16138 




2 


13907 


9488 


-1 


16531 


10060 


8 


14000 


8844 


5 




10233 


11 


12999 


8269 


14 


12730 


10083 


16 


13156 


7410 


19 


12710 


9928 


24 


12687 


6572 


25 


12114 


10085 


28 


12331 


6361 


28 


11874 


10107 


31 


12167 




33 


11066 


10922 


33 


12192 




40 


10747 




42 


11465 




51 


10058 
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Table 5: Measured values of M^ ax , Am 15 , ft(Si II), «i (Si II), v for SN 1999ac and SN 1999aa. 



SN 1999ac 


SN 1999aa 


M max 


-18.98(0.39) a 


-19.14(0.78) 


Amis 


1.30(0.09) a 


0.85(0.08) 


ft(Si II) 


0.098(0.030) 


0.163(0.016) 


«io(Si II) 


8470(300) 


10100(300) 


V 


137(14) 


4(14) 



a As reported in Li et al. (2003). 



